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Целью работы является определение свойств транспортных 
потоков мобильных  строительно-дорожных объектов 
(МСДО) для повышения эффективности управления ими. 
Представлены результаты системного анализа задач 
управления МСДО, изложена концепция их решения и 
выполнена формализация задачи управления движением 
мобильных строительно-дорожных машин в постоянно 
изменяющихся условиях Впервые 


сформулированы системные и структурные критерии оценки 


передвижения. 
состояния мобильных строительно-дорожных объектов, 
проанализированы концепции управления группами 
подвижных объектов, даны рекомендации по выбору типа 
управления. 

транспортных 


Представлены математические модели 
потоков, устанавливающие зависимость 
скорости движения от плотности транспортного потока, а 
также математическая модель оптимизации маршрута 
движения мобильных строительно-дорожных объектов. В 
статье приведены 


интеллектуального 


стратегии адаптивного и 
управления группами подвижных 
объектов различной численности, а также представлены 
результаты регрессионного анализа модели движения 
мобильных объектов. 


Ключевые слова: мобильные строительно-дорожные объек- 
ты, групповое управление, транспортные потоки, организа- 
ция движения, оптимизация маршрута движения. 
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Введение. Вопросы управления подвижными объектами различного вида, безусловно, представляют научный инте- 
рес. При этом особенно остро стоит проблема взаимодействия указанных объектов в условиях недетерминированной 


среды функционирования. Сложившаяся в настоящее время ситуация требует улучшения качества работы транспорт- 
ных комплексов, повышения средней скорости движения объектов, сокращения времени перемещения транспорта 
между пунктами назначения. Решению указанных задач способствует использование информационно-управляющих 


систем. 


Учитывая темпы развития строительной отрасли, а также масштабы застройки современных мегаполисов, необходимо 
разрабатывать новые интеллектуальные системы управления поставками строительных материалов, обеспечивающие: 
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— эффективное использование транспортных и строительных машин; 
— выбор оптимальных маршрутов движения по основным и вспомогательным дорогам. 
В частности, следует: 
— создавать новые модели подвижных объектов и их групп; 
— разрабатывать оригинальные имитационные модели дорожной сети и самонастраивающиеся алгоритмы управле- 
ния, обеспечивающие при необходимости их быструю перестройку. 

Такие алгоритмы способны изменять маршруты движения, перераспределять транспортные потоки, учиты- 

вать скоростной режим на отдельных участках и оптимизировать транспортные затраты при обслуживании строи- 
тельных объектов. 
Формализация задачи управления мобильными строительно-дорожными объектами. Наибольший интерес пред- 
ставляет задача управления мобильными строительно-дорожными объектами (МСДО). Ее решение позволяет добить- 
ся соблюдения графика доставки конструкций и материалов на строительные площадки и сократить транспортные 
расходы. При этом необходимо выполнить комплекс управленческих задач. В первую очередь речь идет об обеспече- 
нии гибкого планирования маршрутов доставки стройматериалов и конструкций. Особенностью управления движени- 
ем МСДО является возможность изменения маршрута в процессе движения транспорта с учетом текущего состояния 
трассы, транспортного запаздывания и средней скорости движения на отдельных участках. 

Проведенный анализ показал, что при управлении МСДО необходимо решить следующие задачи: 

— контроль скорости и маршрута движения МСДО в реальном времени; 

— сбор и обработка информации об отклонении МСДО от графика и маршрута движения; 

— выбор оптимального маршрута движения и последующая его корректировка; 

— контроль состояния технической среды перемещения объектов, определение препятствий и изменение скорости 
перемещения; 

— формирование законов управления движением МСДО и их корректировка с учетом состояния выбранного маршру- 
та движения. 

Концепция управления подвижными объектами и их группами. Применительно к МСДО используются централи- 
зованный или комбинированный принципы управления. Формирование управляющих воздействий для подвижных 
объектов основывается на мониторинге их состояния [1]. Использование различных стратегий управления подвижны- 
ми объектами и их группами позволило установить следующее. В случае применения централизованной стратегии 
при росте численности объектов в группе значительно увеличивается время принятия решений (в сравнении с децен- 
трализованной стратегией, рис. 1). 





Рис. 1. Зависимость времени принятия решения при различных стратегиях управления 


При централизованном управлении решения о действиях объектов группы принимаются одним вычислителем, что 
требует значительного времени на обработку данных (рис. 1, график 1). При этом в соответствии с увеличением числа 
объектов в группе экспоненциально возрастает сложность задачи, а также время ее решения. Децентрализованное 
управление позволяет каждому подвижному объекту принимать самостоятельное решение, обмениваясь информацией 
с другими объектами группы для оптимизации действий. Время управления при росте числа подвижных объектов 
увеличивается постепенно ( рис. 1, график 3). Очевидно, что при смешанной стратегии управления увеличение време- 
ни решения будет демонстрировать некоторое промежуточное значение (рис. 1, график 2). Представленные зависимо- 
сти позволяют обоснованно выбирать тип стратегии управления МСДО. В первую очередь следует определить номи- 


* 
нальное время обслуживания объектов в группе Т., (рис.1). Как видно из рис. 1, при количестве объектов в группе 


меньше № предпочтение следует отдавать централизованной стратегии. Если число объектов в группе превышает М, 
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предпочтение следует отдавать децентрализованной стратегии управления. Смешанную стратегию целесообразно 
использовать при количестве объектов от № до Л. 

Системные и структурные критерии оценки состояния подвижных объектов. Управление МСДО невозможно без 
учета свойств и параметров состояния объектов. Кроме того, необходима формальная постановка задачи группового 
управления [2]. Рассмотрение вопросов группового управления подвижными объектами требует исследования показа- 
телей качества управления и невозможно без учета особенностей технической среды функционирования. Опишем 
состояние подвижных объектов, техносреды и группы объектов с помощью вектор-функции. 

Состояние одиночного подвижного объекта М представим вектор-функцией: 


т 
м =[т(0,т, (0,..т,(0] , (0 
где т; (1) — переменные текущего состояния объекта, определяемые видом объектов. 
Состояние среды, в которой функционирует подвижный объект М, опишем вектор-функцией: 


Е(о = [а (0),е, 0..0], (2) 
где е; (Г) — переменные состояния среды, определяемые видом рассматриваемых задам. 

В качестве переменных состояния среды рассматриваются параметры транспортных потоков, координаты 
препятствий на маршруте, коэффициенты взаимной корреляции плотности потоков на параллельных полосах движе- 
ния, средние скорости движения транспортных потоков и др. 

При групповом управлении МСДО кроме вектор-функций М(!) , описывающих поведение каждого объекта, 


введем описание состояния группы, состоящей из К взаимосвязанных подвижных объектов. Для практических целей 
взаимодействие группы объектов лучше описывать вектором состояния в момент времени [ вида: 


№0) = [1 (0,7, (0..1 0], (3) 
где с; (2) = переменные СОСТОЯНИЯ группы ПОДВИЖНЫХ объектов, определяемые их назначением и решаемыми задача- 
МИ. 

При составлении графика обслуживания строительных объектов определяется плановое время, за которое 
строительная машина проходит весь маршрут Ти и его отдельные участки, совпадающие с участками других маршру- 


тов группы: АТ, (7 = 1,") ь 


где и — число совпадающих участков. 
Основным критерием качества управления группы объектов в данном случае является соблюдение графика 
движения транспортных средств. 
Для полного описания состояния подвижных управляемых объектов необходимо формализовать цель управ- 
ления в виде вектор-функции управляющих воздействий, учитывающих динамику подвижных объектов и техниче- 
ской среды их функционирования: 


т 
О = Фи, и]. (4) 
Используя описанные векторы состояния (2), (3) и (4), состояния одиночного объекта М(Р) представим в ви- 
де: 


м -Е, [соо 0@, С, Е] , (5) 


где С(Р),Е(Р) — динамические составляющие группового поведения объектов и технической среды их функциониро- 


вания. 

Таким образом, состояние подвижных объектов определяется не только управляющими и возмущающими 
воздействиями, но и является результатом динамического взаимодействия рассматриваемых объектов и транспортных 
потоков среды. Под влиянием управления 0(1) объекты могут изменять свое состояние М(1) и состояние группы в це- 
лом. Изменения состояния объектов в общем случае носят непрерывный характер и описываются дифференциальны- 
ми уравнениями: 


№ =Е[МО,Е®, 00] . (6) 
В общем виде задача управления подвижным объектом заключается в определении последовательности его 


состояний М(В, которая обеспечивает достижение сформулированной цели функционирования объекта. Критерием 
поиска оптимальной последовательности действий является экстремум функционала 


У = М, (Еь, №) + | ЕСМ, Е, О@), В , (7) 


1 


| 


где 1, & — начальный и конечный момент времени. 
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Использование данного функционала позволяет учесть влияние среды и оценить качество процесса управле- 
ния. В случае подвижных объектов, выполняющих технологические операции, с помощью сформулированного функ- 
ционала можно определить оптимальный маршрут движения $(1) или оптимальные действия объекта. 
Математическое моделирование маршрутов движения МСДО. Задачи моделирования движения МСДО разделя- 
ются на две основные группы: 

— моделирование плана местности и планирование маршрута движения в соответствии с заданным графиком 
выполнения работ; 
— моделирование динамики движения МСДО при отработке спланированного маршрута. 

Необходимость в таких моделях возникает при создании адаптивных и интеллектуальных систем управления 
строительно-дорожными машинами. 

При моделировании движения важным является отражение в моделях транспортной сети технических средств 
организации движения. Предпочтительным решением этой задачи является разработка специальных модулей 
технических средств и способа их введения в модель транспортно-дорожной среды [3]. При этом технические 


средства организации движения описываются множеством А = {К,} ‚ содержащим в себе подмножества светофорных 
объектов К, ={",}, дорожных знаков К, ={",}, линий дорожной разметки К ={",}. Причем эти объекты 


расположены на дугах е, е Е модели транспортнои среды: 





УкеД(!) 3! е ЕЁ: [оса е,. (8) 
При дислоцировании в модель на участке расположения дуги технического средства этот объект становится 
весовой характеристикой дуги: 
Ут еЁ(!),У е. ЕЕ, т 15 [осшей е <> 39, =®(!: 





(#15 юсшеа ® ‚) л (е, 15 юсшеа ®,) ^( к 15 асйоп е, ), (9) 
где ©(Г) — вектор-функция весовых характеристик. 


Зона действия таких модулей, вводимых в состав модели транспортной среды, ограничивается сечением, 
перекрестком, перегоном, магистралью или всей областью функционирования мобильного объекта. 

Как показали результаты наблюдений за МСДО и результаты моделирования процессов их движения, для 
управления подвижными объектами более адекватной является математическая модель, реализующая зависимость, 
представленную на рис. 2. 





Рис. 2. Взаимозависимость скорости и плотности транспортного потока 


В этом случае зависимость между скоростью и плотностью описывается соотношением [3]: 


—а? —а? 
ехр о а_| (10) 


тах тах 


о=о 


тах 


где о! и 092 — динамические коэффициенты среды, определяемые на основе экспериментальных исследований; 'Оеасы 
Ч их — максимальные значения скорости и плотности потока. 


Тогда интенсивность Движения будет равна: 
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—_ 72 —_ 72 
ехр и р (11) 


=4о 
т й 


тах 


шах тах 


Начальный участок модели адекватен реальному процессу транспортного потока, когда объекты движутся 
свободно и расстояние между ними больше динамического минимума, обеспечивающего безопасность движения и 
необходимую маневренность транспортных единиц. 

Если речь идет о вынужденном движении с синхронизацией скоростей подвижных объектов, целесообразно 
описывать транспортный поток на основе разрывных моделей Ромпома — Блимера: 


4,(а„„-а) 
=.“ 54, <а<ав,, (12) 
Я.А ыы 


тах 





где а, —_ плотность транспортного потока, соответствующая пропускной способности. 


Учитывая разновидности транспортных единиц, в эту модель следует ввести динамические коэффициенты 
ПОДВИЖНыЫхХ объектов. Используя понятие динамической безопасности, введем динамические коэффициенты ДЛИНЫ И 
ширины подвижного объекта: 

| Ь, 


=4 у =-8, 13 
у, 1. Уь ь, (13) 


где /, и, — динамические длина и ширина объекта; /; и ВБ; — статические длина и ширина объекта. 
При этом динамические параметры /, и В, зависят от скорости движения — следовательно, динамические 


коэффициенты длины и ширины также являются переменными и определяются скоростными характеристиками пото- 
ка. Кроме того, значения этих коэффициентов требуют нормализации в зависимости от габаритов мобильного строи- 
тельно-дорожного объекта. 


Оптимизация маршрута движения МСДО. При управлении МСДО с произвольной маршрутизацией одной из 
главных задач является выбор оптимального маршрута с использованием объектных и темпоральных моделей 
территории [4]. Решение задачи должно основываться на задании начального и конечного пункта [5], [6], [7]. Кроме 
того, возможно задание одной или нескольких промежуточных точек по указанию диспетчера. После ввода данных 
алгоритм управления должен определить кратчайший маршрут, отобразить его на карте территории и выдать 
описание маршрута. Для поиска оптимального маршрута лучше использовать ориентированный граф, включающий 
множество альтернативных маршрутов, построенный на основе объектной модели территории (рис. 3) [8]. 


[10 1 Ро [19 Ре 





Ри Раб 


Рис. 3. Ориентированный граф альтернативных маршрутов движения 


Каждую дугу между узлами характеризует вес — длина пути между соответствующими узлами. В соответ- 
ствии с длиной дуги определяются результирующие длины альтернативных маршрутов: 


В о ие (14) 


А>В А>В А>В 


Оптимальным для проезда считается маршрут, на котором обеспечивается пт 5, . Однако при таком подходе 


не всегда обеспечивается минимизация времени обслуживания [9], [10]. Поэтому мы предлагаем после определения 





оптимального маршрута В >В >В >В, >В, —>Р. провести анализ состояния транспортной сети и обстановки 


на выбранном маршруте, выполнить его корректировку с учетом заторов на дороге (рис. 3). Корректировка выполня- 
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ется на основе критерия минимизации времени проезда между пунктами Б, и В ‚ В качестве дополнительной инфор- 


мации используются данные мониторинга состояния дорог. 

Следует отметить, что на работу группы МСДО на маршруте влияет ряд внешних факторов. Для исследования 
этих воздействий используется модель маршрута движения объекта, формируемая на основе модели транспортной 
сети. Для моделирования динамики движения объектов к модели маршрута подключаются источники возмущающих 
воздействий, которые в случайные моменты времени формируют потоки объектов на перекрестках и в начале марш- 
рута. При этом каждому объекту случайным образом задается время появления на перекрестке и полоса движения. 
Кроме того, можно задавать интенсивность потоков на каждой полосе. 

Для оценки характеристик транспортного маршрута сети с нечеткими пропускными способностями следует 
рассматривать ее в виде несвязанного графа, каждая дуга которого характеризуется пропускной способностью 4(хь х)). 
Началом сети является вершина-источник хо = Ро; вершина-сток — х„ = Ре. 

Присутствие нечетких параметров в рассматриваемых задачах обусловлено неопределенностью, 
возникающей под воздействием различных факторов внешней среды (в частности, непредвиденных ситуаций на 
маршрутах). В результате анализа движения МСДО в транспортных потоках можно выделить следующие 
разновидности задач: 

— определение максимального потока транспортной среды с нечеткими пропускными способностями; 
— определение минимальной стоимости маршрута с нечеткими данными; 

— оценка максимального динамического потока с нечеткими данными; 

— задача о потоках с усилениями с нечеткими данными. 

При решении задач определения максимального потока транспортной сети интересны пропускные способно- 

сти, представленные в виде: 
— нечетких чисел [А-типа, 
— интервалов. 
Нечеткие пропускные способности в виде нечетких чисел [К-типа задаются: а, = (0, 4,;0 В,), в, где 


9, +В, =4,. 
Распределение функции принадлежности нечеткого числа имеет вид: 
14 <4 
и, =1 0; 4, $4,4, (15) 
0,4, >а,; 


В задаче о нахождении нечеткого максимального потока речь идет о потоке 4, = {а, Е КОсх,) Е А) от иИСТочЧ- 


ника к стоку. В интервале [0; 1] функция не обязательно линейна и выполняются условия: 


У, Хх, =5 
уеЁ; 4, <а, >(хх,)еА; У а,- 3 Чав=4-ух= (16) 
х,ЕГ(х,) хуЕГ (х,) 0, х. Ех/(5,Ё) 


Первое условие уеЁ означает, что значение у имеет наиболыную степень принадлежности к Ё, где Е 
представляет собой набор значений потока (в пределах интервала), которые допускает эксперт в каждом конкретном 
случае. Второе условие в (16) также означает, что поток 4, меньше либо равен 4, с наибольшей степенью принад- 


лежности. Большее значение функции принадлежности Це (у) показывает, что лицо, принимающее решение, допускает 


большее значение потока у. Значение Е представляет собой нечеткое число [.-типа: Ё = (», +00; ак; 0) ‚в - Выраже- 
Е МЕ 


ния (16) интерпретируются следующим образом: необходимо найти поток, имеющий значение у, причем условие со- 
хранения потока &, < 4, не должно быть нарушено. В работе для решения этой задачи выбран подход, основанный на 


агрегации целевой функции и ограничений. Предложенный алгоритм решения позволяет найти поток, имеющий 
наибольшую величину и одновременно наиболее удовлетворяющий нечеткому ограничению на пропускные способ- 


ности. Нечеткие пропускные способности задаются в виде интервалов: а = [4, а, |, где дь, 4, — соответственно левые 


и правые границы пропускной способности маршрута. 
Алгоритм планирования оптимального маршрута включает следующие основные этапы поиска: 
— преобразование графа на основе контурной модели; 
— нахождение опорных вершин, через которые проходит маршрут объекта; 
— минимизация пути движения с использованием градиентного метода на основе линейных ориентиров. 
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На первом этапе строится матрица смежности А = {а} „», элементы которой представляют собой время про- 
хождения подвижным объектом расстояния от одной вершины графа 1-й области до другой вершины /-й области через 
общую границу по формуле 

а, =5,/у,, (17) 


где а. — вес дуги, п едставляющий собой время п охождения отрезка между вершинами [и |: 5. 
й ы 


у расстояние от 1-и 


вершины до вершины /-й области; у, — скорость движения между 1-й и/-й вершинами. 


Для поиска пути на графе можно воспользоваться модифицированным алгоритмом Дейкстры, в котором ис- 
пользуются индексно-матричный метод поиска и метод градиентной оптимизации маршрута движения объекта. В 
качестве входных данных рассматриваются координаты начала и конца маршрута. В качестве критерия оптимально- 
сти используется минимум времени перемещения. В ходе решения задачи оптимизации определяется маршрут дви- 
жения, состоящий из набора дуг, связывающих определенные вершины графа. 

Целевая функция для оптимизации маршрута движения имеет вид: 


1„„(р,)=У.165;) > ши, (18) 

1 . 
где Р — вершины графа транспортно-дорожной среды, характеризуемые координатами (х,,у,) расположения 
площадей и перекрестков; 5, — дуги графа между вершинами Ри Р; [(5,) — время проезда 5, расстояния от 


(1 1)-й до 1-й границы. 
Время 1,(5,) определяется через расстояние 5 ‚ И среднюю скорость движения потока на данном участке 
Г, : 5,)=5, /У,. 
В качестве ограничений целевой функции следует использовать минимально допустимую скорость движения 
потока на участке при выборе маршрута: 
К, (5,) =У, У, 


19 
В, (5) =У, (5), и 


где у @ Ув (5,) —Щ_ минимальное и максимальное допустимые значения скорости на рассматриваемом участке 


движения 5, ограничивающие 1-ю компоненту целевой функции сверху и снизу. 
Учитывая, что ограничения имеют вид неравенств, воспользуемся методом проекции градиента. В этом слу- 
чае на каждом А-м шаге итерации выполняется процедура возврата компоненты вектора у,(А) на допустимое множе- 


ство значений 5, если у,(К) #5 ‚ т. е. не удовлетворяет условию ограничений (19). Это реализуется за счет присвоения 


у, (К) одного из граничных значений: У» Или У. 


т 


При планировании маршрута движения в качестве начального приближения /„(0) принимается одно из 
направлений смежных дуг рассматриваемой транспортной среды при предельных скоростях движения: 


(0) =Г „(Уи ) ЕМУ 1 (У) Е. (20) 


Градиент целевой функции определяется значениями каждой компоненты вектора, зависит от положения со- 
ответствующей границы и в аналитическом виде представим выражением 


аг„ (5) В Ах’ + Ду’ а А + Ду’, 


ор 


2 2 
45; У, -Д5. Ут" 45°. 








; (21) 


где Ах и Ду — при ащение координат на участке движения объекта. 
м, 


Для реализации пошаговой работы алгоритма необходимо воспользоваться методом градиентной оптимиза- 
ЦИИ Флетчера — Ривса и методом проекции градиента. На каждой итерации работы алгоритма осуществляется про- 
верка выполнения критерия останова, корректировка в случае необходимости шага поиска и проверка на зациклива- 


ние. В качестве критерия останова выбирается неравенство |5 (К) —5(А -1| < &, где  — допустимая ошибка поиска. 
Алгоритм реализует одномерный поиск методом золотого сечения по параметру ^, ‚ минимизирующий функ- 


ЦИЮ 
(2%) =1.„ (5+. 8, ) > шт, (22) 


й 


где &, — антиградиент, начальное приближение которого задается в виде &, = -У/.,„ (5 (0)) : 
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Приближение на шаге итераций вычисляется по формуле 5(&+1) = 5(&)-+^, -Е, ‚ а направление поиска опре- 
деляется соотношением 
[Ул (5%+Б 
[Ми ЗОО | 


Как показали исследования описанного алгоритма, при его работе возможно зацикливание при подходе к точ- 


(23) 





6, =-У1,„ (5+1) +Е, 


ке минимума 5 (К) . Это связано с неправильным выбором № — большое значение способствует зацикливанию про- 


цесса поиска вокруг точки минимума 5 (К) й Для исключения зацикливания программы поиска необходимо ввести 


дополнительную проверку на зацикливание. 
Исследования алгоритма поиска маршрута движения мобильных строительно-дорожных объектов показали 

следующее. Данный алгоритм, исходя из заданного графа транспортной сети, обеспечивает нахождение кратчайшего 

маршрута движения, удовлетворяющего заданным ограничениям по скоростному режиму. 

Выводы. Для подтверждения эффективности использования разработанной модели проведены экспериментальные 

исследования, целью которых было уточнение параметров, характеризующих свойства транспортных потоков МСДО. 

В результате эксперимента получена временная зависимость величины /,(1), которая приведена на рис. 4. 





4 8 12 16 20 24 28 32 56 


Рис. 4. Временная зависимость величины я 5 (1 ) 


Анализ временного ряда подтвердил нормальный закон распределения. Использование вариационного ряда 

позволило получить смешанную модель авторегрессии, описываемую выражением 
К@=0,258+0,789. К а-—П +а()+0,092.а1-—1), (24) 
где а(Г) — белый шум. 

Достоверность модели подтверждена малыми значениями стандартных ошибок параметров модели, которые 
не превышают 6,2 %. Сравнение результата прогнозирования с экспериментальными данными позволило установить, 
что среднее значение доли по данным моделирования составляет 0,2692, а вычисленное на основе экспериментальных 
данных — 0,2703. При этом среднеквадратичное отклонение в первом случае равно 0,1349, а во втором — 0,1213. 

Полученные значения подтверждают высокое соответствие результатов моделирования реальным данным. 
Следовательно, модель авторегрессии может быть использована для прогнозирования условий проезда подвижными 
объектами регулируемых перекрестков на маршруте движения. 
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